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Abstract 

We formulate a new class of stochastic partial differential equations (SPDEs), named 
high- order vector backward SPDEs (B-SPDEs) with jumps, which allow the high-order 
integral-partial differential operators into both drift and diffusion coefficients. Under cer- 
tain type of Lipschitz and linear growth conditions, we develop a method to prove the 
existence and uniqueness of adapted solution to these B-SPDEs with jumps. Comparing 
with the existing discussions on conventional backward stochastic (ordinary) differential 
equations (BSDEs), we need to handle the differentiability of adapted triplet solution 
to the B-SPDEs with jumps, which is a subtle part in justifying our main results due 
to the inconsistency of differential orders on two sides of the B-SPDEs and the partial 
differential operator appeared in the diffusion coefficient. In addition, we also address 
qq ' the issue about the B-SPDEs under certain Markovian random environment and employ 

a B-SPDE with strongly nonlinear partial differential operator in the drift coefficient to 
l/" - ) ' illustrate the usage of our main results in finance. 

o ' 

Key words and phrases: Backward Stochastic Partial Differential Equations with 
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1 Introduction 

Motivated from mean- variance hedging (see, e.g., Dai [10]) and utility based optimal portfolio 
choice (see, e.g., Becherer [3], Musiela and Zariphopoulou [22]) in finance, and multi-channel 
(or multi-valued) image regularization such as color images in computer vision and network 
application (see, e.g., Caselles et al. [6], Tschumperle and Deriche [33, 34, 35], and references 
therein), we formulate a new class of SPDEs, named high-order vector B-SPDEs with jumps, 
which allow high-order integral-partial differential operators C and J into both drift and 
diffusion coefficients as shown in the following equation (1.1), 

r T f T 
(1.1) V{t,x) = H{x) + C{s-,x,V,-)ds+ {j{ S -,x,V,-)-V{s-,x))dW{s) 

t r . 

/ V(s ,x,z,-)N{\ds,x,dz). 
t Jz>0 
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where the operator C depends not only on V, V, V but also on their associated partial deriva- 
tives, i.e., for each integer k > 2 and m > 0, C and J are defined by 

C(s,x,V,-) = C(s,x,V(s,x),V^(s,x),V( S ,x),V^(s,x),V(s,x),-), 

J(s,x,V,-) = (Ms, x,V, ■),..., J d (s,x,V,-)), 

Ji{s,x,V r ) = J l (s,x,V(s,x),V^\s,x),-), ie{l,...,d}. 

Under certain type of Lipschitz and linear growth conditions, we prove the existence and 
uniqueness of adapted triplet solution (V, V, V) to these B-SPDEs. When the partial dif- 
ferential operator C depends only on x, V, V, and V but not on their associated derivatives 
and C = 0, our B-SPDEs with jumps reduce to conventional BSDEs with jumps (see, e.g., 
Becherer [3], Dai [10], Tang and Li [32]). 

BSDEs were first introduced by Bismut [5] and the first result for the existence of an 
adapted solution to a continuous nonlinear BSDE was obtained by Pardoux and Peng [26]. 
Since then, numerous extensions along the line have been conducted, such as, Tang and Li [32] 
get the first adapted solution to a BSDE with Poisson jumps for a fixed terminal time and 
Situ [31] extended the result to the case where the BSDE is with bounded random stopping 
time as its terminal time and non-Lipschitz coefficient. Currently, BSDEs are still an active 
area of research in both theory and applications, see, e.g., Becherer [3], Cohen and Elliott [8], 
Crepey and Matoussi [9], Dai [10], Lepeltier et al. [18], Yin and Mao [36], and references 
therein. 

The study on SPDEs receives a great attention recently (see, e.g., Pardoux [25] and 
Hairer [14]). Particularly, Pardoux and Peng [27] introduces a system of semi- linear parabolic 
SPDE in a backward manner and establish the existence and uniqueness of adapted solution 
to the SPDE under smoothness assumptions on the coefficients, and moreover, the authors in 
[27] also employ backward doubly SDEs (BDSDE) to provide a probabilistic representation 
for the parabolic SDE. Since then, numerous researches have been conducted in terms of weak 
solution and stationary solution to the semi-linear SPDE (see, e.g., Bally and Matoussi [2], 
Zhang and Zhao [38], and references therein). However, our B-SPDEs exhibited in (1.1) 
are fundamentally different from the SPDEs as introduced in Pardoux and Peng [27] and 
as studied in most of the existing researches in the following aspects: First, our system 
formulation is a direct generalization of the conventional BSDEs, i.e., both the drift and 
diffusion coefficients of our B-SPDEs depend on the triplet (V, V, V) and its associated partial 
derivatives not just on V and its associated partial derivatives; Second, our B-SPDEs are 
based on high-order partial derivatives and are subject to jumps. One special case of our 
B-SPDEs available in the literature is the one derived in Musiela and Zariphopoulou [22] for 
the purpose of optimal-utility based portfolio choice, which is strongly nonlinear in the sense 
that is addressed in Lions and Souganidis [19]. 

Note that the B-SPDEs presented in (1.1) are vector B-SPDEs with jumps, which are 
motivated from various aspects such as multi-channel image regularization in computer vision 
and network application through vector PDEs (see, e.g., Caselles et al. [6], Tschumperle 
and Deriche [33, 34, 35], and references therein), coupling and synchronization in random 
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dynamic systems through vector SPDEs (see, e.g., Mueller [21], Chueshov and Schmalfufi [7], 
and references therein). 

To show our formulated system well-posed, we develop a method based on a scheme 
used for conventional BSDEs (see, e.g., Yong and Zhou [37]) to prove the existence and 
uniqueness of adapted solution to our B-SPDEs with jumps in (1.1) under certain Lipschitz 
and linear growth conditions. One fundamental issue we need to handle in the method is the 
differentiability of the triplet solution to our B-SPDEs with jumps, which is a subtle part 
in the analysis due to the inconsistency of differential orders on two sides of the B-SPDEs 
and the partial differential operators appeared in the diffusion coefficient. So more involved 
functional spaces and techniques are required. In addition, although there is no perfect 
theory in dealing with the strongly nonlinear SPDEs (see, e.g., Pardoux [25]), our discussions 
about the adapted solution to (1.1) can provide some reasonable interpretation concerning 
the unique existence of adapted solution before a random bankruptcy time to the strongly 
nonlinear B-SPDE derived in Musiela and Zariphopoulou [22]. 

In the paper, we also provide some discussion concerning our B-SPDEs under random 
environment, e.g., the variable x in (1.1) is replaced by a continuous Markovian process X{-). 
To be convenient for readers, we present a rough graph in Figure 1 with respect to sample 
surfaces for a solution to a B-SPDE and in terms of sample curves for a solution to the 
B-SPDE under random environment. 




Figure 1: Sample surfaces and sample paths for a B-SPDE and a B-SPDE under random environment 

The rest of the paper is organized as follows. In Section 2, we first introduce a class of 
B-SPDEs with jumps in finite space domain, then we state and prove our main theorem. 
In Section 3, we extend our discussions in the previous section to the case corresponding to 
infinite space domain and under random environment. In Section 4, we use an example to 
illustrate the usage of our main results in finance. 
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2 A Class of B-SPDEs with Jumps in Finite Space Domain 



2.1 Required Probability and Functional Spaces 

First of all, we introduce some notations to be used in the paper. Let (0, J 7 , P) be a fixed 
complete probability space on which are defined a standard d-dimensional Brownian motion 
W = {W(t),t G [0,T]} with W{t) = (Wi(t), W d (t))' and h- dimensional subordinator 
L = {L(t),t G [0, T]} with L(t) = (Li(t), Lft(t))' and cadlag sample paths for some 
fixed T G [0, oo) (see, e.g., Applebaum [1], Bertoin [4], and Sato [30] for more details about 
subordinators and Levy processes), where the prime denotes the corresponding transpose of 
a matrix or a vector. Moreover, W, L and their components are assumed to be independent 
of each other. In addition, each subordinator Lj with % G {l,...,h} can be represented by 
(see, e.g., Theorem 13.4 and Corollary 13.7 in Kallenberg [17]) 

(2.1) Li(t) =Oit+ f I ZiNi(ds,dzi), t > 

J(0,t] Jzi>0 

where Ni((0,t] x A) = J2o<s<t Ia(L(s) — Lj(s - )) denotes a Poisson random measure with 
a deterministic, time-homogeneous intensity measure dsvi{dzi), where Ia(') is the index 
function over the set A, the constant ai is taken to be zero, and V{ is the Levy measure. 
Related to the probability space (£2, J 7 , P), we suppose that there is a filtration {J~t}t>o 
with F t = o{W(s),L{\s) : < s < t} for each t G [0,T], A = (Ai,...A fc )' > 0, and 
L(Xs) = (Li(Aia),...,L fc (A h a))'. 

Secondly, let AT = {1, 2, } and D be a close connected domain in BP for a given p G TV. 
Then we can use C k (D, R q ) for each k,p,q £ Af to denote the Banach space of all functions 
/ having continuous derivatives up to the order k with the following uniform norm, 

(2 - 2) mck ^ = c^imwmiS W C)(X) 

for each / G C k (D,R q ), where r(c) for each c G {0,1,..., fc} is the total number of the 
following partial derivatives of the order c 

(c) 9 C f r (x) 

(2 ' 3) = dx n[ dx h 

with G {0, 1, c}, Z G {1, r G {1, q}, and «i + ... + i p = c. Moreover, for the late 

purpose, we let 

( 2 - 4 ) f{h,...,i P ) = (fuh,...,i P )'-"'fq%i,...,i P )^ 

(2-5) fM(x) = (fi°\x),... y f% ) (x)), 

where each j G {l,...,r(c)} corresponds to a p-tuple (ii,...,i p ) and a r G {1, ...,(/}. In 
addition, let C°°(D,R q ) denote the following Banach space, i.e., 

(2.6) C°°{D,R?) = {fe flC fc ( J D, J R"),||/|| c .o 0(Dig) <oo| 
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where 



(2.7) 



k=l 



for some discrete function with respect to k G {0, 1,2, ...}, which is fast decaying in k. For 
convenience, we take £(&) = e~ k . 

Thirdly, we introduce some measurable spaces to be used in the sequel. Let Lj-QO, T]; R q ) 
denote the set of all i? 9 - valued measurable stochastic processes Z(t,x) adapted to {J~t,t G 
[0, T]} for each x G D, which are in C°°(D, R q ) for each fixed t G [0, T]), such that 



(2.8) 



E 



l^(*)llc ,00 (D, ? )^ 



< OO 



and let Ljr p ([0, T], R q ) denote the corresponding set of predictable processes (see, e.g., Def- 
inition 5.2 and Definition 1.1 respectively in pages 21 and 45 of Ikeda and Watanabe [15]). 
Moreover, let Ljr T (Q;R q ) denote the set of all i? 9 -valued, .Fr-measurable random variables 
£(x) for each x G D, where £(x) G C OD (D,R q ) satisfies 



(2.9) 



E 



< OO. 



In addition, let L^QO, T], R h ) be the set of all i? fe -valued predictable processes V(t,x,z) 
(Vi(t,x,z), Vh{t,x,z))' for each x G D and z G R+, satisfying 



(2.10) 
and let 



E 



h r j< 

E 



,i=i 



Jzi>0 



Vi(t-,z) 



C°°(D,q) 



Vi(dzi)dt 



< OO 



(2.11) Ll c (DxR h + ,R q * h ) = \v:DxR h + ^R q * h ,J2 f H^llc^,,) < oo} 

I i=i J ^>o J 

with the associated norm for any v G E 2 UC {D x R+, R qxh ) and c G {0, 1, oo} as follows, 

\\v\kc= J z> Jvi(zi)\\ 2 cc{Dq) \m{dzi^ . 



(2.12) 



In the end, we define 

(2.13) Q 2 A[0,T]) = L 2 T ([0,T],R q ) x L% p ([0,T], R qd ) x L 2 p ([0, T], R qxh ). 

2.2 The B-SPDEs 

First of all, we introduce a class of g-dimensional B-SPDEs with jumps and terminal random 
variable H(x) G Ljr T (Q;R q ) for each x G D as presented in (1.1), where for each s G [0, T] 
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and z = (zi, ...,z h ) G R\, 

V(s,-) = (V 1 (s r ),...,V d (s,-))€C°°(D,W* d ), 
V(s,-,z) = (V 1 (s,;z 1 ),...,V h (s,;Z h ))eC 00 (D,R^ h ), 
N (Xds , x , dz) = (Ni(Xids,x,dzi), Nh(Xhds,x,dzh)Y ■ 

Moreover, in (1.1), C is a (/-dimensional integral-partial differential operator satisfying, a.s., 

(2.14) A£( c \s,x,u,v) < K D (\\u - v\\ ck+c{DA) + ||fZ - v\\c™+°{p,qd) + II s ~ v\\v,c) 

for any (u, u, u), (v, v, v) G C k (D, R?) x C m (D, Ri xd ) x L 2 uc {R h + , Ri xh ) with c G {0, 1, oo}, 
where Kjj depending on the domain D is a nonnegative constant, \\A\\ is the largest absolute 
value of entries (or components) of the given matrix (or vector) A, and 

(2.15) A£( c \s ,x,u,v) = £ {c \s,x,u,-)-C^{s , X, v, ■). 

Similarly, J = (J~i, ...,Jd) is a q x (i-dimensional partial differential operator satisfying, a.s., 

(2.16) \\AJ^(s,x,u,v)\\ <K D {\\ 
Moreover, we suppose that 

(2.17) C ic) (s,x,u, •) < K D (||u||c*+c(£>,,) + \\v\\c™+°(D,qd) + l|w|L,c) ' 

(2.18) J^(s,x,u,-) < K D \\ 

u \\c m + c (D,q)- 

Example 2.1 The following conventional linear partial differential operators satisfy the con- 
ditions as stated in (2.14.)- (2.18), 



(Cu)(t, x) 



d 2 u{t,x 
dxjdxj 



<9u(i, a;) 



+ c{x)u(t, x) 



+ X M 3 



d 2 u(t, x) 



chx(i, a;) 



+ c(x)u(t, x) 



(Ju)(t,x) 



J2 M x ) 



d 2 u(t,: 



d 



dx j dx q 

1,3=1 J 3=1 



du(t, x) 
dxj 



+ c(x)u(t, x), 



where a\j (x), bj (x), c^ c \x), a\j (x), bp{x), and S- C \x) are uniformly bounded over all 
x G D and i,j G {1, d} and c G {0, 1, 2, ....}. 

Theorem 2.1 Under conditions of (2.14)-(2.18), if £(t,x,v,-) and J(t,x,v,-) are {Tt}- 
adapted for each fixed x G D and any given (v, v, v) G C°°(D, R q ) x C°°(D, R qxd ) x L 2 u oo (D x 
R\,Ri xh ) with 



(2.19) 



£(-,x,0,-),J(-,x,0,-) GL^([0,T],^), 
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then the B-SPDE (1.1) has a unique adapted solution satisfying, for each x G D and z G R\, 
(2.20) (V(-, x), V(-,x), V(.,x, z)) G Q%([0, T}) 

where V is a cadlag process and the uniqueness is in the sense: if there exists another solution 
(U(t,x),U(t,x),U(t,x,z)) as required, we have 



E 



f T (\Mt) - V{t)\\c-(D, q ) + Wit) - V(t)\\ 2 c~ { D, qd ) + W{t) - n*)ll*,oc) dt 
J o 



0. 



We divide the proof of Theorem 2.1 into the following three lemmas. 

Lemma 2.1 Under the conditions of Theorem 2.1, for each fixed x G D, z G R+, and a 
triplet 



(2.21) 



(U(,x),U(;x),U(,x,z)) G Q 2 A[0,T}), 



there exists another triplet (V(-, x), V(-, x), V(-, x, z) such that 

(2.22) V(t,x) = H(x)+ C(s-,x,U,-)ds+ {J(s~,x, U,-)- V(s~, x)) dW(s 



t Jz>0 



V(s ,x,z)N(\ds,x,dz), 



where V is a {Ft}- adapted cadlag process, V and V are the corresponding predictable pro- 
cesses, and for each x G D, 



(2.23) 
(2.24) 
(2.25) 



E 



E 



E 



\\V{t,x)\\ 2 dt 



T 



\\V{t,x)\\ 2 dt 



< CO, 

< oo, 



E 

i=l ' 
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/7 


Vi(t ,x,z) 


Vi(dzi)dt 


< oo 


/0 Jzi>0 









Proof. First of all, for each fixed x G D, z G R+, and a triplet (U(-,x), U(-,x), U(-,x,z)) 
as stated in (2.21), it follows from conditions (2.14)-(2.19) that 

(2.26) C(;x,U,-)eL 2 T ([0,T],R«), J(, x, U, •) G L 2 T ([0, T], ^ xd ). 

Now consider C and J in (2.26) as two new starting C(-,x, 0, •) and J{-,x, 0, •), then it 
follows from the Martingale representation theorem (see, e.g., Lemma 2.3 in Tang and Li [32]) 
that there exists a unique pair of predictable processes (V(-,x), V(-,x, z)) which are square- 
integrable for each x G D in the senses of (2.24)-(2.25) such that 



(2.27) V(t,x) = E 



H(x)+ f T ' C(s~,x,U,-)ds+ [ T J(s-,x,U,-)dW( 
Jo Jo 



= V(0,x) + f t V(s-,x)dW(s) + f [ V{s~,x,z)N{Xds,x,dz) 

Jo Jo Jz>0 
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which implies that 



(2.28) V(0,x) = H(x)+ / C(s-,x,U,-)ds+ / J(s~, x, U, -)dW(s 



V(s~,x)dW(s) - 



V(s ,x,z)N(\ds,x,dz). 
o Jz>0 



Moreover, due to the Corollary in page 8 of Protter [29], V(-,x) can be taken as a cadlag 
process. Now we define a process V as follows, 



(2.29) V(t,x) = E 



H(x)+ [ T C(s-,x,U,-)ds+ [ T J(s-,x,U,-)dW(s) 
Jt Jt 



Then by simple calculation, we know that V(-,x) is square-integrable in the sense of (2.23), 
and moreover, it follows from (2.27)-(2.29) that 



(2.30) 



V{t,x) = V(t,x) - /* £(s-,x,U,-)ds - f J(s-,x,U,-)dW(s) 
Jo Jo 



which indicates that V(-,x) is a cadlag process. Furthermore, for a given triplet (U(-,x), 
U(-,x), U(-,x), U(-,x,z)), it follows from (2.27)-(2.28) and (2.30) that the corresponding 
triplet (V(-,x), V(-,x), V(-,x,z)) satisfies the equation (2.22) as stated in the lemma, which 
also implies that 

(2.31) V(t,x) = V(0,x) - f f jC{s-,x,U,-)ds - [* (J(s~,x,U,-) - V (s~ , x)) dW (s) 



o 



+ 



V(s ,x,z)N(Xds,x,dz). 



10 Jz>0 

Hence we complete the proof of Lemma 2.1. □ 



Lemma 2.2 Under the conditions of Theorem 2.1, for each fixed x G D, z G R+, and a 
triplet as in (2.21), we define V(t,x), V(t,x), V(t,x,z) through (2.22). Then (V^X^x), 
V(°\-,x), V( c \-,x, z)) for each c G {0, 1, } exists a. s. and satisfies a. s. 

(2-32) Vg^frx) = H^ t Jx) + J^^ tp) ( S -,x,U,-)ds 



[ T [ V® .Js-,x,z)N(Xd8,x,dz), 

Jt J z>0 



i.[S 



where i\ + ... + i v = c and i\ G {0,1,..., c} with I G {l,...,p}. Moreover, i for each 



c G {0,1,...} is a {Ft^-adapted cadlag process, 



(ii-.-ip) 



and Vh 



(c) 



(h—ip) 



(ii—ip) 
are the corresponding 



predictable processes, which are square-integrable in the senses of (2.23)- (2.25). 
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Proof. First of all, we show that the claim in the lemma is true for c = 1. To do so, for 
each given t G [0, T], x G D, z G and (U (t, x), U (t, x), U(t, x, z) as in the lemma, let 



(2.33) 



(V^\t,x),vgHt,x),vg\t,x,z)) 



be defined through (2.22) where C and J are replaced by their first-order partial derivatives 
£y) and jQ) in terms of Xi with I G {l,...,p}. Then we can show that the triplet defined 
in (2.33) for each I is indeed the required first-order partial derivative of (V, V, V) that is 
defined through (2.22) for the given (U,U,U). 

As a matter of fact, for each / G {U, U, U, V, V, V, N}, small enough positive constant 5, 
and I G {l,...,p}, define 

(2.34) f {l):S (t,x) = f(t,x + S ei ), 

where e; is the unit vector whose Ith component is one and others are zero. Moreover, let 
(2.35) 



A.). ^",- /ft ' l -A.) 



W v ' ,5 ■'(0 

for each / G {17, 17, L>, V, V, V}. In addition, let 

(2.36) AZ^ s (s,x,U) = ^(l(s,x + 8e l ,U(s,x + 8e l ),-)-l(s,x,U(s,x),-)) 

-lj^(s,x,U(s,x),-) 

for each 1 G {£, J}. Then, by applying the Ito's formula (see, e.g., Theorem 1.14 and 
Theorem 1.16 in pages 6-9 of 0ksendal and Sulem [24]) to the function 



aAV ( ^Jt,x)) = Tr[AV ( )tUt,x))e 



(i),s 



(0,<5 V 



for some 7 > 0, where Ti(A) denotes the trace of the matrix A' A for a given matrix A, we 
have 

(2.37) C(AyW(t,,))+ [ T Tr(AJ^ s (s,x,U)-AV$/ Sl x))e^ds 

J t 

+ It L >o ^ ( A ^(S (S ~' X ' Z) ) e2rS ti(^ds, x, dz) 
= 2 jT (-7TY (AF ( « (s,x)) + (Ay«( S ,x))' (A£^( S , x, [/))) e 2 ^s - M(i) 

^ Tr (A^(s, s)) e 27S ds + 7 jf T || A^jJ^a, x, C/)|| 2 e 27S ds - Af (t) 



= 7 / T A£g 5 ( S , x, 17) 2 e^ds - Mi(t) 

*/ 1 

if, in the last equality, we take 
(2.38) 



~ 37^ 2 
7 = -r^ > 0, 



2 7 
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where M§(t) is a martingale of the following form, 

d ,. T 
It 



2 E r (^ s (s-,x))' (A(40jg,( fl - *,C0 - A(^)« 5 ( S -,x)) e^W^s) 



j=l ,Jt Jzj>0 N 7 x 7 

+ 2 E / T / (^i^"^))' ((VW)(i),«(*~ e^ 8 {Nj) {l) ,s{\jds,x,dzj). 

Now, it follows from Lemma 1.3 in pages 6-7 of Peskir and Shiryaev [28] that, for each 
t G [0, T] and a > 0, there is a sequence of {<5 n , n = 1,2, ...} C [0, a] such that 



(2.39) 



= E 



ess sup C( A vl 1 ^ (<,*)) 

0<<5<cr V h 



ess sup a*V$ Jt,x)) 

{5„:0<<5„<<7,n=l,2,...} V 7 



lim £7 

n— »oo 



C(Ay ( «jt,x)) 



< 7 lim I? 

. _ n 

/•T 

< jE ' 



Atfl Sn (s,x,U) 



lux L E[M Sn (t)} 



ess sup 

0<<5<ct 



e 2 ^ s ds 



where "esssup" denotes the essential supremum and the first inequality in (2.39) follows from 
(2.37). So, by the Lebesgue's dominated convergence theorem, we have 



(2.40) 



lim E 



ess sup C(AV$ s (t,x)) 



0<8<a 



< jE 



lim ess sup A£$ x (s, x, U) 
t <r->o o<6<* {l) ' S 



' e 2 ^ s ds 



since, due to the mean- value theorem and the conditions stated in (2.17), we have 



A£W s (t,x,U) < 2K D ( \\U\\ C H + i (D>q) + \\U\\ C m +HD , qd) + 



U 



Then, by (2.40) and the Fatou's lemma, we know that, for any sequence a n satisfying a n — > 
along n G AT, there is a subsequence TV 7 C AT such that 



(2.41) 



ess sup C(AV/jj ^(t, x)) — > along n £ M' a.s., 



0<<5<<t„ 



which implies that the first-order derivative of V in terms of x\ for each / G {1, exists 
and equals V^\t,x) a.s. for each t G [0,T] and x G D, and moreover, it is {J" t }-adapted. 
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Next, it follows from the similar proof as used in (2.39) that 

rT 



(2.42) 



lim E 



' ess sup Tr (AjJJ^, x, U) - AV$ s (s, xfj e 2ls ds 



< lE 



[ T lh 



lim ess sup 

It °"~>0 0<<5<cr 



AC 



(i) 

(l),5 



[s,x,U) 



^ s ds 



Thus, by (2.41) and (2.42), we know that 



lim Av£\(t, x) = lim Aj£\(t, x,U) = a.s. 

which implies that the first-order derivative of V with respect to x; for each Z G {l,...,p} 
exists and equals V"i^(i,x) a.s. for every £ G [0, T] and x £ D, and moreover, it is a {^"t}- 

predictable process. Similarly, we can get the conclusion for V^\t, x, z) associated with each 

/ 2 1 * X - Z • 

Secondly, assuming that (y( c_1 )(i, x), F( c_1 )(i, x), y( c_1 )(i, x, 2)) corresponding to a 

given (U(t,x), U(t,x), U(t,x,z)) G QjKt ) 7 !) exists for an y S iven c G I 1 ; 2 ,---}- Then 
we can show that 



(2.43) 



(y( c )(t,x),y( c )(t,x),y( c )(t,x,z)) 



exists for the given c G {1,2, ...}. 

As a matter of fact, consider any fixed nonnegative integer numbers ii,...,i p satisfying 
ii + ... + i p = c — 1 for the given c G {1, 2, ...}, each / € {V, V, V}, each Z G {1, ...,p}, and 
each small enough <5 > 0, let 



(2.44) 



(c-l) 



correspond to 

1 i p )( s ' x + 5eu C/ ( s ' x + 6ei ^ ') with X G { jC > ^ via ( 2 - 22 )' where the de- 
ferential operators £ and J" are replaced by their (c — l)th-order partial derivatives £^~ 1 ] ) 

and 4i"t)- Similarl y' let ( F (S..(i i +i)...i p )(*' a; )' ^(ii ) ...(<«+i)-<p)f t ' x )' ^(S..(i J +i)...ip)(*' a ;' 2! )) be 

defined through (2.22) where C and J are replaced by their cth-order partial derivatives 
£it-(ii+i)-i p and ^h...{n+i)...i p corresponding to a given t,x,U(t,x), U(t,x), U(t,x,z). More- 
over, define 



(2.45) 



A f(h...(i l +i)...i p ),s( t > x ) ~ 
for each / G {U, U, U, V, V, V} and let 

1 

6 



fjl.X+V-ip)/^^ ^ifj _ -(c) /. x 



= ^( I (t 1 t)( t ' a; + (5e '^( t ' a; + (5e ')'-) 

X £T.t) ( s > x > ^ x )> •)) " I £i...(i J +i)...i p ) (*, f/(s, x), •) 
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for X G {£, J}. Then, by the Ito's formula and repeating the procedure as used in the second 
step, we know that (vfi r , ^ • \(t,x), vfi ,. . \(t,x), vfi /■,-,>, ■ \(t,x,z)) exist for 
the given c G {1, 2, ...} and all / G {1, .. .,£>}, which implies that the claim in (2.43) is true. 

Thirdly, it follows from the induction method with respect tocG{l,2,...} that the claims 
stated in the lemma are true. Hence we finish the proof of Lemma 2.2. □ 

Lemma 2.3 Under the conditions of Theorem 2.1, all the claims in the theorem are true. 

Proof. Let D 2 r ([0,T],R q ) be the set of .Revalued {Jf}-adapted and square integrable 
cadlag processes as in (2.8). Moreover, for any given 7 G R, define M®{0, T] to be the 
following Banach space (see, e.g., the similar explanation as used in Yong and Zhou [37], and 
Situ [31]) 

(2.47) M°[0,T\ = D 2 T ([0,T],Ri) x Lj^ p ([0, T],R qxd ) x L 2 p ([0,T], R q * h ) 

endowed with the norm: for any given (U, U, U) G M®{0, T], 



(2.48) 



(U, U, U) 



M 



D = X)e(fc) (U,U,U) 
k=l 



7,fc 



where, without loss of generality, we assume that m = k in (1.1) and 



(2.49) (U,U,U) 



M D . 

1,K 



E 



sup \\U(t)\\ 2 ck(Dg) e 

0<t<T 



27* 



+ E 



f T \\u(t)\ 

Jo 



C k {D,qd) 



^dt 



+E 



ii{t) 



^dt 



In addition, through (2.22), we can define the following map, 

H : (U(-,x),U(;x),U(;X,z)) -> (V (■ , x) , V (• , x) , V (• , x, z)) . 

Then, based on the norm defined in (2.48), we can show that S forms a contraction mapping 
in M%[0,T\. As a matter of fact, consider (£/*(-, x), tP(; x), IP(; x, z)) G M%[0,T\ and 
(V^(-,x), ^(-,x,z)) = £(£/*(-, x), 17*(-, x), [/*(•, x,*)) with i G {1,2,...}, define 

Afi = f i+i_ f i with /G |^^ ^ y ^^j 

and take 

(2.50) C(ACT(t, x)) = Tr (AET(t, x)) e 27 '. 

Then it follows from (2.14) and the similar argument as used in proving (2.37) that, for a 
7 > and each i G {2, 3, ...}, 

- 1 )-A^(s,x))e 27S ds 



(2.51) 



C(AE/*(t,s)) + 

rT 



Tr(AJ(s,x, 
Tr(AL^(s _ ,x,z)) e 27S A>(Ads,a;,(iz) 
< 7^' ||A£(s,x, U\ U' 1 - 1 )^ e 2 ~< s ds - M (i) (t) 
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where K a is some nonnegative constant depending only on Kp, and for the last inequality, 
we have taken 



(2.52) 



3K 



1 



D 



2 7 



> 0. 



Moreover, N l (t) appeared in (2.51) is given by 



(2.53) N^it) 



r 



ACT" 1 (a) 



C k (D,q) 



+ 



ACT" 1 to 



C k {D,qd) 



+ 



AW'^s) 



v.k 



e 2 ^ s ds 



and M l (t) is a martingale of the following form, 



(2.54) M\t) 

" 2 E f T ((A^)( S -,x))' (AJjis-^lP,^- 1 ) - (AC/^Or, s)) e 2 ~< s dW 3 {s) 
j=i Jt 

+2E^ T ^ >o ((AC/^( S -,x))'((C/ i ) j ( S -, a ;,^)) e^'NjiXjds^dzj). 

j = l 3 

Then, by (2.51)-(2.54) and the martingale properties related to stochastic integral, we have 



(2.55) 



E 



AU\t,x)\ e 27 * + y Tr (^A J (s, x, U\ W' 1 ) - AW (s,x))e 2 ^ds 



AW(s) 



e 2 ~< s ds 



<j(T + l)K a (AW' 1 , AW' 1 , AW' 1 ' 



7,fc 



Next it follows from (2.54), the Burkholder-Davis-Gundy's inequality (see, e.g., Theorem 
48 in page 193 of Protter [29]) that 



(2.56) E 



sup 

0<t<T 



M\t) 



< 

3=1 



sup 

0<t<T 



J q (AW(s~,x))' 



+4j> 

3=1 
d 

3=1 
h 

+K b Y / E 

3=1 



sup 

0<t<T 



f I (Af/ i (s-,x)) / (AC/ i ) i (s-,x,z i )e 27S iV(A j (is,x,(iz j 

JO Jzi>0 V ' 



AW(s,x) (AJ^^x.UW*- 1 ) - {AW) 3 {s,x) 



e^ s ds 



in 

2 



AW(s,x) (ACT)j(s,x,^) e^'XjVjfadzrfds) 



2 
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< K b E 



sup \\AU\t,x)\\ z e 

v O<t<T 



2 2 7 t 



< -£ 
~ 2 



sup ||A^(t,x)|| 2 e 27 ' 

0<<<T 



+dK 2 E 



+K 2 E 



Jo 



^ s ds 



T 



AU\s) 



j/,0 



* 2* 



sup \\AU\t,x)\\ 2 c0{q) e 2 ^ 

0<t<T 



+ j(T + l)dK^E N l ~ l (t) 



where K b is some nonnegative constant depending only on and T, and we have used 
(2.55) for the last inequality of (2.56). Thus it follows from (2.51)-(2.56) that 



(2.57) 



E 



sup 

0<t<T 



AU\t) 



,2~ft 



C°(q) 



< 2 (l + dKl) 7(T+1) \{AU' 1 - 1 , AU' 1 - 1 , AU 1 ' 1 ) 



M D , 



Moreover it follows from (2.51) and (2.17) that, for i G {3,4, ...}, 
(2.58) E Tr ^A^7*(s, x)j e 2ls ds 

< 2E J\r (A^s,!,^,^ 1 ) - AU\s,xj) e 2 ^ s ds 
+2E J^Tr [AJ{s,x,U i ,U i - 1 )) e 2 ^ s d 

< 2^K C ( (ACT -1 , ALT -1 , ALT -1 



7,fc 



+ 



(AC^ -2 , ALT" 2 , ALT -2 ; 



7,fe . 



where iT c is some nonnegative constant depending only on Kp and T. Thus it follows from 
(2.51) and (2.57)-(2.58) that 

(2.59) \\(AU\AU\AW)f MD ^ 



(Atr- 2 ,Atr- 2 ,Atr- 2 ; 



M 1 ? 

7>& / 



14 



where is some nonnegative constant depending only on and T. 

Now, by Lemma 2.2 and the similar construction as in (2.50), for each c £ {1,2,...}, we 
can define 



(2.60) 
where 



C{AU c '\t,x)) = Tr (aC/ c '^,x)) e 27 ', 



AU°>%x)) = (AU(°)' < (t,s)),Al7W' < (*,x)),...,Al7( c )' < (t,x)) / . 
Then it follows from the Ito's formula and the similar discussion for (2.59) that 

(2.61) (AC/*, AC/*, ACT*) ' 



7<c 



< 7 K d [ (ACT -1 , ACT -1 , ACT -1, 



which implies that 

(2.62) (AW, AW, AW 



7,fc + c 



+ 



(ACT -2 , ACT -2 , ALT -2 ; 



< 7^ ( (AIT*- 1 , ACT" 1 , ALT" 1 ) 



+ 



(Atr-^A^-^A^- 2 ) 



Since (a 2 + 6 2 ) 1 / 2 < a + b for a, 6 > 0, we have 



(2.63) 



(Air, ait, ait) 



< ^^[(ACT^AC/*- 1 ^*- 1 )^^ 



where Kf is some nonnegative constant depending only on Kn and T. Therefore, by taking 
7 > large enough such that 2y /A fKf sufficiently small and by (2.62), we know that 



(2.64) 



]T II (ACT*, ACT*, ACT* 



i=3 



< 



< oo. 



(ACT 2 , AU 2 , ACT 2 



+ 



(ACT 1 , ACT 1 , ACT 1 ) 



Thus, from (2.64), we see that (U t ,U' l ,U l ) with i G {1,2,...} forms a Cauchy sequence in 
A4®{0, T], which implies that there is some (CT, U, U) such that 

(2.65) (CT*,CT*,CT*) (U,U,U) as i -»• oo in M°[0,T]. 



Finally, by (2.65) and the similar procedure as used for Theorem 5.2.1 in pages 68-71 of 
0ksendal [23], we can finish the proof of Lemma 2.3. □ 

Proof of Theorem 2.1. 

By combining Lemma 2.1-Lemma 2.3, we can reach a proof for Theorem 2.1. □ 
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3 B-SPDEs in Infinite Space Domain and under Random Environment 
3.1 B-SPDEs in the Infinite Space Domain 

First of all, for a given nonnegative integer 6 and each n G {6+1, 6+2, ...}, define the following 
sequence of sets 

(3.1) D n = {x G R p : b < \\x\\ < n} 
and let 

(3.2) R p = {x G R p : \\x\\ > 6} 

Moreover, let C°°(R p ,q) be the Banach space endowed with the following norm 

oo 

(3-3) WfWc^RU) = E Z(n + l)\\f\\c-(D n , q) 

n=b+l 

for each / G C°°(Rl, q), and let 

(3.4) Q 2 A[0,T}) = L 2 A[0,T],R«) x L% p ([0,T],RO d ) x L 2 p ([0,T], R«* h ) 

be the corresponding space defined in (2.13) when the norm in (2.7) is replaced by the 
associated one given in (3.3). 

Theorem 3.1 Assuming that there exists a nonnegative constant K R p such that conditions 
(2.14)-(2.18) are satisfied when Kd is replaced by K R v . Moreover, if £(t,x,v, •) and^f(t,x,v,-) 
are {Tt}- adapted for each x G R p , z G R\, and any given (v,v,v) G C°°(R P , R q ) x 
C°°(R p b ,R qxd ) x Ll j00 (R p x R^,Ri xh ) with 

(3.5) £(-, x, 0, •), J(;x, 0, •) G 1% ([0, T],R q ) , 
then the B-SPDE (1.1) has a unique andapted solution satisfying, 

(3.6) (V(;X),V(;X),V(;X,Z)) G 

where V is a cddldg process and the uniqueness is in the sense: if there exists another solution 
(U(t,x),U(t,x),U(t,x, z)) as required, we have 



(3.7) E 



T 



(\\U{t) - V{t)f Coo{Rlq) + \\U(t) - V{t)f Coo{RUd) + \\U(t) - V(t)\\l,oo) dt 



0. 



Proof. It follows from (3.5) and the similar argument used for (2.62) in the proof of 
Lemma 2.3 that 

(3.8) (U\;X),U\;X),U\;X,Z)) G Q^([0,T]) 
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with (U°,U°,U°) = (0,0,0), where (U 1 , U 1 , U 1 ) is defined through (2.22) in Lemma 2.1. 
Then, over each {D n } with n £ {b + 1,6 + 2,...}, it follows from (2.62) in the proof of 
Lemma 2.3 that 



(3.9) 



{ait, ait, Air 



Mr 



<7K g [ (ACT -1 , AC/ 1 " 1 , AC/ 



H— 1 Ar~r»— 1" 



Mr 



(AC/* -2 , AC/' -2 , AC/ 1 



where -KT 9 is some nonnegative constant depending only on T and B%. Then it follows from 
(3.9) that the remaining proof for Theorem 3.1 can be conducted similarly as in the proof of 
Theorem 2.1. Hence we finish the proof of Theorem 3.1. □ 

3.2 B-SPDEs under Random Environment 

Assuming that the random environment under consideration is characterized by a i? p -valued 
Markov process X(-) with continuous sample paths and its associated stopping time r is 
defined by 



(3.10) 



T = mf{t G [0,T], \\X(t)\\ < b}. 



Then the g-dimensional B-SPDEs with jumps under random environment X(-) can be de- 
scribed as follows, 



(3.11) 



V(t,X(t)) 



H(X(t)) + £ C(s-,X(t),V,-)d8 
+ jT (j(s-,X(t),V, •) - V(s-,X(t)) dW{s) 
- T f V{s-,X{t),z,-)N{\ds,X{t),dz). 

Jt Jz>0 

Moreover, define 

(3.12) &A[0,t]) = ^([0,r]^) x L^([0,r],^) x L 2 ([0, r], i?^) 

be the corresponding space defined in (2.13) when the norm in (2.7) is replaced by the 
following one, 



(3.13) 



||/(x)||oo = sup /^(a 



where f^\x) for each x € R? is defined as in (2.5). 

Theorem 3.2 Under the conditions as stated in Theorem 3.1, the B-SPDE in (3.11) under 
random environment X(-) has a unique adapted solution satisfying, 



(3.14) 



{V{.,X{.)),V{;X{.)),V{;X{-),z)) G (&([0,r]), 
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where V is a cadlag process and the uniqueness is in the sense: if there exists another solution 
(U(t,X(t)),U(t,X(t)),U(t,X(t))) as required, we have 



E 



+ 



(\\U(t,X(t)) - V(t,X(t))\\l + \\U(t,X(t)) - V(t,X(t))\\ 
\\U(t,X(t))-V(t,X(t))\\l :00 )dt 



2 

oo 



Proof. First of all, it follows from the similar discussion as in Situ [31] that Qf([0, t]) is a 
Banach space. Then we know that all of the claims stated in Theorem 3.1 are true over the 
space Qf([0,t]), which imply that the claims in the current theorem are true. □ 

Example 3.1 The solution V(t,x, •) to the B-SPDE in (1.1) is described by random surfaces 
and the solution V(t, X(t)) to the B-SPDE under random environment in (3.11) is represented 
by random paths, which are shown in Figure 1 presented in the Introduction. 

4 An Illustrative Example in Finance 

In this section, we consider a financial market consisting of two assets and an external random 
factor. One asset is supposed to be a risk-free account whose price So(t) is subject to the 
following ordinary differential equation, 



(4.1) 



dS (t) = rS (t)dt, 
S o (0) = 1, 



where the interest rate r is a nonnegative constant. Another asset is stock whose price process 
S(t) satisfies the following SDE for each t G [0,T], 



(4.2) 



dS(t) = S(t)(3(Y(t))dt + S{t)a(Y(t))dWi{t), 
S(0) =s>0, 



where the random factor Y(t) with t G [0,T] satisfies 
(4.3) 



dY(t) = c(Y(t))dt + d(Y(t)) (pdW^t) + y/T=p*dW2(tj) , 
y(0) = y€R 



with p G (—1, 1). Moreover, we suppose that the market coefficients f = f3,a,c,d satisfy the 
standard global Lipschitz and linear growth conditions and a(y) > k > for all y G R and 
some positive constant k. 

Beginning at t = with an initial endowment x G R+ , an investor invests at any time t > 
in the risky and riskless asset. The present value of the amounts invested are denoted, respec- 
tively, by 7i"o(t) and vri(t), and then the present value of the investor's aggregate investment 
is given by X w (t) = 7To(£) + ni(t), which satisfies (see, e.g., Musiela and Zairphopoulou [22]) 

(4.4) dX w (t) = cr(t)7r(t) • (\(t)dt + dW{t)) 



18 



where " • " denotes the inner product, Tr(t) = (iro(t), tti dW = (dW\, dW2)' , and 
(4.5) m -(W®\- T ^' 



\ *(»(<)) 

Moreover, for a given constant b > 1, let the following r be the bankruptcy time for the 
investor, 

(4.6) r = inf{t > 0,X n (t) < b}. 

One objective to study the above financial system is to find the optimal portfolio choice 
based on maximal expected utility of terminal wealth over all admissible strategies (see, e.g., 
Merton [20]), i.e., to solve the following stochastic dynamic optimization problem, 

(4.7) V(t,x) = sup E p [u t (X*(t))\ F t ,X"(t) = x] 

Ar 

where A T denotes the set of all admissible strategies tt: ir(t) is self-financing and {Ft}- 
progressively measurable, satisfying 



E 



(ir(s)a(s)) 2 ds 



< oo and X n (t) > with t € [0,r], 



and the utility is taken to be the following constant relative risk aversion (CRRA) case: 

u T (x) = — with < 7 < 1, 7 / 0. 

7 

Then it follows from the discussions in Musiela and Zariphopoulou [22], 0ksendal [23], 0ksendal 
and Sulem [24] that the value function process defined in (4.7) should satisfy the following 
B-SPDE, 

(4.8) V M = M X(r))+ f^WW^ 

Jt ^v xx {s,x) 

+ J\v(s,x))'dW(s) 

where a(t) = (a(Y(t)),0)', {a{t)) + = (l/a{Y(t)),0), and dW = (dW 1 ,dW 2 y. 

As pointed out in Musiela and Zariphopoulou [22], the B-SPDE in (4.8) is newly derived 
and belongs to a class of strongly nonlinear B-SPDEs (see, e.g., the related discussion in Lions 
and Souganidis [19]). However, based on Theorem 3.1 in the previous section of the current 
paper and the discussion in Musiela and Zariphopoulou [22], we can show that there exists a 
unique adapted solution, before a random bankruptcy time (i.e., over [0,r]), to the B-SPDE 
in (4.8) over the class of functions satisfying the conditions required by Theorem 3.1. In fact, 
based on the discussions in Musiela and Zariphopoulou [22], Glosten et al. [13], we have the 
following observation that there is a pair of V and V satisfying (4.8), i.e., 

(4.9) V(t,x) = -x'f(t,Y(t)) 5 

7 
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where / is a solution of the following partial differential equation 



ft + 7; d (y)fyy + (civ) + — — 




i>?{y)f 

25(l- 7 ) 



= 



with f(r,y) 



1 and 5 = (1 — 7)/(l — 7 + p 2 7), and 



(4.11) 



(4.10) 



Vi(t,x) 



V 2 (t,x) 




Thus it follows from (4.9)-(4.11) that (V, V) is a solution to (4.8) such that the conditions 
in Theorem 3.1 are satisfied over [0, r] and hence it is the unique adapted solution to (4.8) 
over [0,t]. Moreover, it follows from Theorem 3.2 in the previous section and the discussion 
in Musiela and Zariphopoulou [22] that, for each t < s < r, the optimal feedback portfolio 
process is given as follows, 
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